Backaction limits on self-sustained optomechanical oscillations 
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The maximum amplitude of mechanical oscillators coupled to optical cavities are studied both 
analytically and numerically. The optical backaction on the resonator enables self-sustained os- 
cillations whose limit cycle is set by the dynamic range of the cavity. The maximum attainable 
amplitude and the phonon generation quantum efficiency of the backaction process are studied for 
both unresolved and resolved cavities. Quantum efficiencies far exceeding one are found in the re- 
solved sideband regime where the amplitude is low. On the other hand the maximum amplitude is 
found in the unresolved system. Finally, the role of mechanical nonlinearities is addressed. 

PACS numbers: 85.85.+j, 42.50.Wk 



I. INTRODUCTION 

Optomechanics is a rapidly growing field in which the 
interaction between optical cavities and mechanical res- 
onators is studied, see Refs. [fi]-|} for recent reviews. A 
lot of interest lies in fundamental questions such as the 
role of gravity in quantum decoherence, the ultimate lim- 
its on position detection, and backaction evasion. These 
questions might be answered by cooling the mechanical 
resonator to its ground state. Only very recently this has 
been achieved using the photon pressure in microwave 
and optical cavities |j) . Using red detuned light the ther- 
mal vibrations of the resonator could be cooled below the 
zero-point motion, which typically lies in the femtomcter 
range Q. On the other hand, for practical applications, 
one would like to have amplitudes as large as possible. 
This can be achieved with blue detuned light, which am- 
plifies the motion of the resonator. For example, we have 
recently demonstrated a non- volatile mechanical mem- 
ory j?J and synchronization of remote mechanical oscil- 
lators H using the large motion amplitudes generated 
by the optical backaction of an on-chip racetrack cav- 
ity. Large amplitude self-sustained oscillations (SSOs) 
also enabled the observation of chaotic dynamics H and 
the zero- frequency anomaly . Reaching high amplitude 
motion is thus important for both technological advances 
as for fundamental research. Here, we address the ques- 
tion what ultimately limits the maximum amplitude of 
regenerative oscillations in a cavity-optomechanical sys- 
tem. 

Figure |l|a shows a schematic of a Fabry-Perot cavity 
where one of the mirrors can move, forming an optome- 
chanical resonator. The analysis presented here is not 
limited to this particular system, but can be applied 
to any cavity-optomechanical system including racetrack 
cavities [Q, photonic crystal structures pj| , and 
even to microwave-cavity optomechanical systems fl2| - 
All these systems have in common that a displace- 



ment changes the cavity frequency and thereby the num- 
ber of photons inside it. This results in a backaction 
on the mechanical element. In the following, first the 
cavity-resonator dynamics will be introduced in Sec. ||, 
followed by a study of the appearance of self-sustained 
oscillations (Sec. III). Section IV shows what happens in 
the large amplitude regime when the harmonic approxi- 
mation breaks down. Finally, Sees. [v| and VI study the 
quantum efficiency and the role of mechanical nonlinear- 
ities respectively. 



II. CAVITY-RESONATOR DYNAMICS 

The threshold for the onset of optomcchanically- 
induced SSO is well studied, both experimentally and 
theoretically [§, ||, ^5[-p0|: For small harmonic motion, 
the cavity field is only slightly perturbed and the photon 
number oscillates at the mechanical frequency with an 
in-phase and a quadrature part as illustrated in Fig. [l]b. 
In this linear regime the cavity response is proportional 
to the motion amplitude. The quadrature part of the 
oscillating photon occupation changes the damping rate 
from its intrinsic value 70: For red detuned light the op- 
tical backaction increases the total damping rate 7, lead- 
ing to cooling |2l|-|2^]. For blue detuning the damping 
is reduced and can even become negative, the so-called 
dynamic instability. In that case the amplitude grows 
exponentially until it becomes limited by nonlinearities 
in either the resonator or in the cavity, and the cavity 
dynamics is strongly perturbed by the oscillator motion 
(Fig. |l|c). The main focus of this work is on the role of 
the cavity, but the question of how mechanical nonlin- 
earities affect the maximum amplitude will be addressed 
in Sec. 



VI 



The coupled equations of motion for the displacement 
u of a harmonic oscillator and the optical field inside 
the cavity \/hflr,a (in the frame rotating at the laser 
frequency, are [[TtJ [l9[ B2l : 
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FIG. 1: (Color online) (a) Schematic representation of a 
generic optomechanical system. Laser light with power P 
enters the Fabry-Perot cavity through the left mirror and 
bounces back and forth between the two mirrors. The right 
mirror acts as a mechanical resonator, (b) Optical backaction 
in the small amplitude limit. A delay between the displace- 
ment u and the change in the photon occupation n due to the 
finite cavity lifetime k _1 leads to ellipsoidal trajectories whose 
area is proportional to the work done on the resonator. Nega- 
tive (positive) work leads to cooling (heating) of the resonator 
in the red (blue) detuned region, (c) The optical backaction in 
the large amplitude regime where the oscillator sweeps across 
the entire cavity resonance. The black line in (b,c) indicates 
the position and linewidth of the cavity, (d) Simulated ringup 
of a mechanical resonator with 70 = 0.001, coupled to a cavity 
with k = 100, Ao = 30, and com = 2000. From this timetrace 
the amplification rate A~ 1 dA/di and oscillation frequency ujr 
are derived (e). 



Here, ujq and m are the resonator frequency and mass, 
and n = a* a is the cavity photon number. The detun- 
ing is the frequency difference between the laser and the 
cavity. For zero displacement this is: Ao = ~ ^c|«=o; 
a displacement changes this to A = Ao + gowiu, where 
goM = —dfl c /du is the optomechanical coupling con- 
stant. When the detuning is zero, the laser with power P 
fills the cavity with photons until the occupation reaches 
the steady-state value n max = 4Pk c / K 2 hflL, where the 
cavity linewidth k is the sum of the external and intrin- 
sic linewidths (k c and ft, respectively). Note that for 
most situations the linear coupling between the cavity 
frequency and the displacement suffices, but refinements 
have been proposed. These include a displacement- 
dependent goM jR, quadratic coupling [fzif , multiple opti- 
cal resonances [£5|, |26) , and effects due to moving bound- 
aries [§7|, |§. 

The full model (Eqs. [j] and ||) contains 8 parameters, 
which are not all independent. By writing the model in 



a dimcnsionlcss form 

u = —u — ^fou + coMCia* (3) 
a = — i(Ao + u)a — ±ka + (4) 

the number of independent parameters is reduced to 4 
p9| . Here, frequencies have been normalized by Uo, and 
displacements by the lengthscale ojo/^om, so that t = uj^t 
and u = goMu/uiQ. Furthermore, t^ 1 = uio/lo is the in- 
trinsic quality factor, and the cavity field is scaled as 
a = a I -y/n max and n = n/n max . The coupling strength 
is c m = 2n max u% pm g£ )M /u$, indicating that the cou- 
pling can be viewed as the ratio of the zero-point-motion- 
induced fluctuations of the cavity frequency compared to 
the resonator frequency. Also note that the optomechan- 
ical coupling coefficient goM appears in com squared due 
to the forward- and backaction 0. 

Equation || can be rewritten in terms of the complex 
amplitude U = (u — iu) exp(— it) [ j30[ by discarding fast 
oscillating terms at frequencies ~ 2cl>o (the rotating wave 
approximation) : 

U = --70U - icoAin^, where h Uo = {h(t)e~ lt }. (5) 

This describes the evolution of the slowly varying am- 
plitude A = \U\ and phase 9 = IU. n Ua is the Fourier 
component of the radiation pressure at the oscillation 
frequency, indicating that only frequencies near the res- 
onance frequency contribute. In the absence of coupling 
to the cavity (i.e., com — 0) A decays exponentially 
back to zero at a rate 570. With coupling present, the 
out-of-phase part — com (n sin(i + 9)} of the photon pres- 
sure changes the amplitude, whereas the in-phase part 
—coM(ncos(t + 9)} changes the phase of the oscillations. 
As explained above, the change in photon number is 
proportional to a small displacement with some delay. 
This means that the cavity response can be written as 
n^n = where the complex response function $ de- 
pends on Ao and h. The real part of <f> is responsible 
for the optical spring effect |32j as it changes 6 at 
a constant rate. On the other hand, the imaginary part 
modifies the damping rate from 70 to 7 = 70 — 2coMlm$. 

III. SELF-SUSTAINED OSCILLATIONS 

From the discussion in the previous section, it follows 
that when the imaginary part of $ exceeds 7o/2com the 
total damping becomes negative and oscillations start to 
grow exponentially as illustrated in Fig. [j]d. Ultimately, 
the oscillations are limited in amplitude since the power 
provided by the cavity is finite and the dynamics reaches 
a limit cycle [^9|. This can be understood as follows: 
when the oscillations become too large, is no longer 
linear in U since the cavity occupation is a nonlinear func- 
tion of u, a Lorentzian to be precise. Thus, when the os- 
cillations exceed the range where only the first order term 
of a Taylor expansion of n around u = is needed (i.e. it 
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exceeds the dynamic range), the proportionality between 
h u and U no longer holds, or, cquivalently, $ becomes 
amplitude dependent. An experimental signature of this 
is the appearance of dips in the optical power coming 
out of the cavity 0, [H| |33|-|3l| which results when 
the oscillations sweep past the cavity resonance peak. In 
principle $ could depend on the complex amplitude at 
all past times which would make the analysis of the sys- 
tem challenging. However, since the optical cavity field 
adapts much faster (~ than the resonator ampli- 

tude changes (~ la 1 ), $ only depends on the present 
amplitude. Moreover, it is independent of the phase of 
the oscillations, 9. $ is thus only a function of A. Fig. 
|l|e shows the evolution of the amplification rate and the 
oscillation frequency, which are directly related to the 
imaginary and real part of $, respectively. After the 
transient (Fig. |l|d and e) the oscillator reaches a steady 
state. The backaction-limited amplitude of this limit cy- 
cle is a solution to A = which is, according to Eq. ||, 
identical to Im$(A) = 7 /2coa/- Note that because the 
absolute phase 9 is not fixed, thermal or backaction force 
noise will lead to a slow phase diffusion of the oscillator 

The conclusion of the above discussion is that, to find 
the maximum amplitude, one needs to know the func- 
tion $(A). When the backaction is not too strong and 
the amplitude is not too high, the motion is to a good 
approximation harmonic. However, as will be shown in 
Sec. IV for large amplitudes this is no longer valid and 
the coupled resonator-cavity dynamics will have to be 
calculated by numerical integration. The cavity response 
for u(t) = Acos(t) was calculated analytically by Mar- 
quardt et al. p9| as an infinite sum of Bessel functions: 
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so that the amplitude-dependent cavity response be- 
comes: 
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(7) 

Figure shows the amplitude-dependence of the mag- 



nitude and phase of $ 



J A for various values of R 



ranging from the deeply resolved to the very unresolved 
sideband regimes. In both regimes, a flat region exists 
on the left side of the plot where the oscillations of n are 
much smaller than n max and the cavity response is linear. 
By taking the limit A —> of Eq. ^ the cavity response 
at this plateau is obtained: 



onto a single curve for the USR. Fig. || also shows that 
the constant region extends up to A ~ R in the USR; for 
larger amplitudes the cavity response rolls off smoothly 
as A' 3 . 

In the resolved sideband regime (RSR) the situation 
is different: still there is a constant region at low am- 
plitudes, but the magnitude of k<& now scales as R 2 
which also follows from Eq. ||. This can be under- 
stood as follows: the cavity linewidth acts as a second- 
order low pass filter for n(t) that filters out fast oscil- 
lations. The slope of 20 dB/dcc means that A ~ 1 
is needed to have a response in n ~ 1. This value of 
A thus demarcates the end of the plateau as shown in 
the Figure. For larger amplitudes, $ is, unlike in the 
USR, not a smooth function of A. It contains many dips 
that become deeper for smaller R as shown in Fig. 

lb- 
Here, only a few terms contribute to the sum in Eq. [?] 
which makes the oscillations in the asymptotic form of 
Jm(A) -> (2/ 7 ri) 1 / 2 cos(i + ?7i7r/2 + tt/4) for A — > oo ap- 
parent. In contrast, in the USR many terms contribute 
and the oscillations arc washed out, resulting in smooth 
curves. 

Since the amplitude of the limit cycle is determined 
by the condition that the total damping rate is zero 
(Im$ = 7o/2coM; see above) the oscillations of &(A) in 
the USR lead to a multitude of solutions for the limit cy- 
cle of the oscillator for resolved optomechanical systems 
p9[ , whereas there is a unique solution for an unresolved 
system. When expanding the equation of motion for A 
(which is obtained from Eq. || for small excursions 8 A 
around the limit cycle with amplitude A) one finds: 
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8 A, j A = -2c m Im 



dl 



(9) 



A=A 



For positive Im[<9$/9A] the limit cycle is unstable and 
the excursions grow until reaching another solution where 
both the damping rate is zero and where, at the same 
time, the derivative is negative. Note that even though 
the limit cycle itself is characterized by a vanishing damp- 
ing rate j(A) = 0, perturbations of the oscillator ampli- 
tude are overdamped and return to A at a rate 7i/2. At 
the fixed point (the limit cycle in the u-ii plane corre- 
sponds to a fixed point in the U representation [p[) the 
oscillation period is ujq(1 — coMRe < ! ) ), which differs from 
that of the uncoupled harmonic oscillator due the pres- 
ence of the nonlinear optical potential p9| . 



IV. THE LARGE AMPLITUDE REGIME 



$(i -> 0) 



A 



[Ag + m 2 + hi*? 



A 2 



(8) 



In the unresolved sideband regime (USR) for large R (and 
a constant ratio between R and Ao), $(A — > 0) oc R^ 1 . 
This is the reason that the plateaus in Fig. collapsed 



In the discussion above it was assumed that the mo- 
tion of the oscillator is harmonic during the entire period. 
This is a good approximation in the RSR, where the cav- 
ity cannot respond quickly enough to the fast mechani- 
cal oscillations and where the amplitudes are relatively 
small. 
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FIG. 2: (Color online) Magnitude (a,b) and phase (c,d) of 
the steady state response of the cavity photon number 
to a harmonic displacement with varying amplitude A for 
Ao + (u) — 0.30ft. The curves with different values for k are 
extracted from numerical time-domain simulations (App. |X|). 
The right panels show zooms of the area around A = 10. 



On the other hand, in the USR the amplitude will be 
shown to be much larger and the cavity detuning os- 
cillates back and forth between large negative and pos- 
itive values. There the photon number thus changes 
rapidly from zero to a large number and then back to 
zero. This kicks the resonator every time the detuning 
crosses zero. When the resonator is moving forward (i.e. 
to the right in Fig. [l]a) work is done on the resonator 
and the gained energy is transferred back to the cavity 
during the backward motion. This is also reflected in 
the dynamics of the oscillator. Instead of the ellipsoidal 
shape for harmonic motion, Fig. ^a shows a mushroom- 
like phase portrait with two different amplitudes {A\ and 
A^). The step in the velocity coincides with large peaks 
in n. Note that both these effects were observed in a op- 
tomechanical system consisting of a Bosc-Einstein con- 
densate inside a Fabry-Perot cavity j3§]. If the oscilla- 
tion would sweep through the resonance slowly, the cavity 
field would always be in equilibrium with the input field 
and the peaks would have the same height n max since 
n(t) = n max /[{2A(u(t))/K} 2 + 1]. The work done is then 
fogoM J_ oc n(t)v(t) dt = = Fj 7m max^ K j where the minus 
and plus signs are for left and right moving resonators, 




FIG. 3: (Color online) Limit cycles of large amplitude motion 
in the USR for 70 = 0.01, R = 100, A = 30 and com = 
2 ■ 10 s . (a) Displacement- velocity phase portrait with the two 
amplitudes (A\ and A2) and the direction of motion indicated, 
(b) Three dimensional phase portrait with the photon number 
on the vertical axis. The black line shows the projection onto 
the u-u plane, (c) Time trace of the cavity occupation. 



respectively. Note that the limits of the integral have 
been extended to ±00 since long before and long after 
the resonance is hit, the cavity detuning is so large that 
the occupation is almost zero at those times. Naturally, 
the work is proportional to the laser power via n max and 
to the linewidth: a larger k makes the Lorentzian line- 
shape wider and hence its area (to which the work is 
proportional) larger. The work is, however, independent 
of the optomechanical coupling coefficient since the <?om- 
term in the force Tigounit) cancels the one originating 
from the detuning A = Ao + <?omm- More importantly, 
since the contributions for the left and right moving tra- 
jectories are equal, the net work during a whole cycle is 
zero. Similar to the sideband formalism for small mo- 
tion |l7], |2j| a delay between the displacement and 
the cavity response is thus also needed in the large am- 
plitude regime to have the net energy transfer required 
to overcome the mechanical damping. One effect of such 
a delayed cavity response is shown in ||b and c. Inter- 
estingly, the timetrace of n{t) shows two different peak 
heights. When the amplitude is small {A\) the velocity 
is also low and the cavity can fill up for a longer time 
than during the backward motion with a larger ampli- 
tude (A2) and a corresponding larger velocity. Finally, 
the fast oscillations at the tail of the peaks in Fig. gc 
result from interference between the input field and the 
Doppler-shifted cavity photons that entered the cavity at 
earlier times [^9|, j35fl . 

The amplitude of the limit cycle is thus determined 
by the balance between the net energy gained during 
one period and the intrinsic damping. Figure ^a shows 
the maximum amplitude for different coupling strengths 
plotted against the cavity linewidth. For low coupling 
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FIG. 4: (Color online) Amplitude (a) and quantum efficiency 
(b) for 70 = 0.01 and Ao = 3, obtained by direct numerical 
integration of Eqs. |^ and ^ (see App. ^) The dotted lines 
in (a) indicate the amplitude predicted by Eq. [ll] and the 
shaded region indicates where A > k and A < k 2 (see text). 



the range over which the SSO are stable is narrow, but 
this increases with increasing com- In the RSR the am- 
plitude is rather small but increases with increasing k. 
The timetraces from which Fig. [| is derived (for details 
see App. |a|) show that in the RSR for large coupling 
( c om <L 100) the system is chaotic (|]. The scatter of 
the curves in this regime are reminiscent of this. Also 
note the multistability of the amplitude in the curve for 
com = 10 |p9[ . The amplitude keeps on increasing with 
« while approaching the USR with an exponent of about 
2/3. Then far in the USR the maximum occurs and is 
finally followed by sharp drop on the right. This divides 
the region with stable SSO from the one where the out- 
of-phase part of the backaction is too small to create a 
negative damping rate. 

From the simulations it is clear that the maximum am- 
plitude occurs in the fiat region in the unresolved side- 
band regime. There, the dynamics corresponds to that 
of Fig. M and the first step to analyze the maximum at- 
tainable amplitude is to consider the asymmetry in the 
peaks. For a given initial amplitude, say A\, the work 
done during the kick can be calculated and from that 
the final amplitude Ai can be obtained by balancing the 
difference in work to the mechanical damping [ ^6| . How- 
ever, even though the peaks have different heights due to 
the velocity difference, the net work done in one cycle is 
still zero when accelerations during the kick are not taken 
into account. To proceed, first it is assumed that the kick 
happens fast so that the displacements before and after 
are close, which is the case when the time spent within 
the cavity lincwidth is small, i.e. Ai ^ S> K. Secondly, it 
is assumed that the detuning does not change too fast so 
that A/k 3 , A 2 /ft 4 and higher order derivatives are small. 
This corresponds to Ai 2 <C ft 2 - The shaded area in Fig. 
Ha shows that both conditions are satisfied for the flat 



region near the maximum amplitude for the larger cou- 
plings. The cavity response then becomes: 



n(A,A) 



1 



1 + (2A/ft) 2 



32AA/k 



[1 + (2A/S) 2 ]' 



(10) 



which has the usual Lorentzian shape for adiabatically 
slow motion that does not yield a net energy transfer, but 
the second term provides a correction for finite velocity. 

When multiplying Eq. [Io| by A and integrating over time 
for both the forward and backward motion, the net en- 
ergy gain is found. Then by setting this equal to the dissi- 
pated energy in one cycle the amplitude A 2 = \ {A\ + A\) 
becomes: 
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(11) 



The dotted lines in Fig. |]a show that the maximum am- 
plitude found from the numerical simulations is indeed 
given by Eq. [ll] when both assumptions are satisfied. For 
smaller coupling strength the shaded area is narrow and 
a deviation between the simulated and calculated am- 
plitude develops, indicating that the assumptions break 
down. However, the simulations show that the amplitude 
is larger than predicted using Eq. [ll], which thus still 
provides a lower bound. Interestingly, Eq. [ll] shows that 
the amplitude of the oscillations is independent of the 
cavity decay rate. Although the work done during each 
kick is oc ft the net energy transfered does not depend on 
k since the latter requires a delay in the cavity response, 
which is cx 1/ft. These two contributions thus balance 
each other. Also A depends only weakly on the mechan- 
ical damping rate 70 and on the coupling strength com- 
Yet the largest amplitudes are obtained for the strongest 
coupling and highest mechanical quality factors. Also, 
although one needs blue detuning to start the SSOs, the 
final amplitude is independent of the detuning. This can 
be understood as follows: when the resonator amplitude 
is so large that it rapidly sweeps over the cavity and the 
decay rate is fast enough to empty the cavity before the 
next kick (c.f. Fig. ||) the exact time at which the kicks 
happen (as set by the condition u{t) = — Ao/<?om) is not 
really important. Finally, note that the 30 dB/dec slope 
of $ shown in Fig. ||a is reflected in the exponent of 1 /3 
in Eq. |TT]. 



V. QUANTUM EFFICIENCY 

In the previous section it was shown how to reach 
the largest amplitudes, however this does not necessarily 
mean that the process is efficient. The quantum effi- 
ciency QE quantifies how many phonons are generated 
by a single photon. Semiclassically, it is the ratio of the 
number of phonons dissipated by the intrinsic damping 
7o x ±itlujqA 2 /huio and the rate of photons entering the 
cavity r„ . Note that the total rate of photons generated 
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in the laser P/MIl is not used since not all of these pho- 
tons enter the cavity: this depends on the ratio of Kj and 
k c . Since the quantum efficiency is an intrinsic property 
of the cavity-resonator system, only the photons that ac- 
tually enter the cavity are taken into account. Using 
Eq. ky the rate equation for the photon number is found: 

1/2 

h = —nn + r in , with r in (t) = Kfj„; ia (a + a*)/2. The time- 
dependence of Tin enters via the non-trivial dynamics of 
a. However, by averaging over one period of the me- 
chanical resonator and noting that for periodic motion 
(h) = 0, the relation (rj„) = nn between the average 
number of phonons in the cavity n and the average rate 
is obtained. The quantum efficiency then becomes: 



nn \2u zpm 



To 



A 2 



(12) 



where the term in the brackets is identified as the number 
of phonons in the mechanical resonator. This shows that 
the larger the amplitude becomes with a smaller average 
number of photons, the larger QE. 

From the simulated timetraces both the motion ampli- 
tude and the average photon number are obtained. Fig- 
ure [pD shows the dependence of the simulated quantum 
efficiency for the data shown in panel (a) on the linewidth 
and coupling strength. In the deeply USR QE is low be- 
cause there the cavity can fill up fast and n is relatively 
large. The large amplitude reached in this regime thus 
requires a large amount of photons and the process is not 
efficient. When decreasing k, A drops (Fig. [|a), but yet 
the quantum efficiency goes up as n drops faster than A 2 . 
Finally, QE saturates to a coM-dependent value above 1 
when approaching the RSR. 

At first it might seem surprising that the quantum ef- 
ficiency in the RSR can exceed one. Note that this is not 
an artifact of the definition of rj n : the same holds when 
using the total laser power. The quantum efficiency is 
most easily understood using the sideband picture of the 
optical backaction Q. Assume that the laser is blue de- 
tuned from the cavity by ~ ujq. Photons cannot enter the 
cavity since their energy is not within the linewidth of the 
cavity. However, by emitting a phonon the photons end 
up in the Stokes (m = —1) sideband and have the right 
energy to enter the cavity. This seems to imply that the 
QE can be at most 1 in the RSR. To understand why 
QE can be much larger than unity, the sidebands (i.e. 
Fourier coefficients) a m of a(t) are analyzed by expand- 
ing the exp(— iAsint) term in Eq. Q resulting in: 



~a{t)e- imi At = - 



2 ^ 

n— — oo 



Jn(A)J n — m (A) 

in + iAo + \k 



(13) 



resonance. In the USR more values of n lie within the 
cavity linewidth and both the m = — 1 and m = +1 side- 
bands are involved in the dynamics (Fig. ||b) . The small 
asymmetry between the Stokes and anti-Stokes sidebands 
provides the energy to the resonator needed to overcome 
the intrinsic damping. Also note the difference in ampli- 
tudes: In the RSR the carrier (m = 0) has a small am- 
plitude as it is detuned from the cavity and the m = — 1 
sideband is of the same order as the carrier. In the USR 
the situation is reversed: Now the laser light is within the 
cavity linewidth and hence the carrier is of order unity, 
but the sidebands are much smaller. This indicates that 
the phonon emission efficiency is lower in the USR com- 
pared to the RSR. 

For small amplitudes there are thus only one or two 
sidebands, but for large amplitudes the situation is dif- 
ferent (Fig. ^c,d): now many more Bessel functions have 
a nonvanishing value and more sidebands emerge. The 
number of sidebands can be estimated as follows: the mo- 
tion of the oscillator frequency modulates the light inside 
the cavity up and down. Light that entered the cavity at 
time to becomes frequency shifted by goM[u(t) — it(to)] 
at time t. The number of available sidebands in the RSR 
is thus roughly equal to the dimensionless amplitude A. 
Figure [|c and d show that is indeed the case. In the USR 
a m decays faster than in the RSR due to the shorter life- 
time. Finally, the overall quantum efficiency is the bal- 
ance between the number of available sidebands and the 
probability to emit a single phonon. The sidebands for 
small motion in Fig. ||a and b showed that the latter 
process is very efficient in the RSR but less efficient in 
the USR, although QE can still exceed one in the USR 
for large coupling. Combining all of this indicates that 
in the RSR the quantum efficiency could be as high as 
QE ~4>1, which is the result of multi-phonon emis- 
sion. This maximum value is indeed reached just before 
the end of the SSOs on the left side of Fig. [|a and b. 



VI. OPTICAL VS MECHANICAL 
NONLINEARITIES 

So far only the optical nonlinearities that limit the mo- 
tion of the self-sustained oscillations have been addressed. 
Another source of nonlinearity that can limit the ampli- 
tude is mechanical in origin: when the motion is large, 
stress induced in the resonator modifies the resonance 
frequency from the small displacement value. For mod- 
erate amplitudes, the equation of motion becomes that 
of the well-known Duffing resonator, now combined with 
optomechanical backaction: 



This shows that when Ao ~ 1 in the RSR a pole oc- 
curs for n = — 1. Now for small A only Jo and J±\ have 
an appreciable amplitude, and hence only the m = — 1 
sideband becomes occupied as illustrated in Fig. |^a. 
This sideband corresponds to the emission of one phonon 
whereby the blue-detuned photon ends up at the cavity 



-70U — loqU — WqCW 3 + hgoMn/m. (14) 



The nonlinearity parameter a has the units of m~ 2 and 
the scaled version is defined as a — «Wq/iJq M . A posi- 
tive value of a indicates a stiffening mechanical spring, 
whereas a negative value corresponds to a spring constant 
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FIG. 5: (Color online) Motion-induced sidebands of the opti- 
cal field for small (top) and large amplitude (bottom) in the 
RSR (left) and USR (right). Negative m corresponds to the 
emission of phonons, positive m to absorbtion. For Ao = +1 
the cavity resonance is located at m — — 1 (dotted lines). 



that decreases with increasing amplitude. The most im- 
portant effect of the cubic term in this context is that it 
modifies the frequency of the oscillations when the am- 
plitude changes. For small a, the oscillator simply oscil- 
lates at the optical-spring-modified resonance frequency. 
However, when a is large the oscillation frequency shifts. 
Figures ^a and b show the oscillation frequency of two 
nonlinear resonators oscillating due to the optical back- 
action of an unresolved (a) and a resolved optical cavity 
(b). In the unresolved case, the oscillation frequency is 
shifted by almost a factor of two for a — 0.002 compared 
to the harmonic oscillator. For negative a the frequency 
can be lowered to about 90% of the original value before 
the resonator escapes from the metastable state around 
u = (sec inset). Fig. ||c shows the dependence of 
the motion amplitude on a. Although there is a depen- 
dence, only a ~ 10% change is visible. This is because 
in the USR the backaction is broadband in nature, and 
a change in resonance frequency will hardly affect the 
amplitude of the oscillator. In the RSR the situation is 
different. Here the laser should be close to the blue side- 
band of the cavity (i.e., at Ao + jom{u) = wr) to have 
strong backaction heating. A change in the oscillation 
frequency with amplitude will move the detuning away 
from this optimal value and thus reduce the resulting os- 
cillation amplitude. The simulations in Fig. [|i show 
that this is indeed the case. The oscillation amplitude 
is strongly peaked around a = and drops when mov- 
ing away from this value. For large magnitudes of a the 
feedback mechanism mentioned above clamps the ampli- 
tude to a value of A ~ 1. Note that this can be used 
to create a tunable optomechanical oscillator. First the 
nonlinear resonator is set into oscillation with the opti- 
mal Ao, but then the detuning is slowly adjusted. Since 
the amplitude is strongly peaked for u>r = Aq + ,9om (u) 



FIG. 6: (Color online) Oscillation frequency (a,b) and am- 
plitude (c,d) of nonlinear resonators in the USR (k = 100, 
left) and RSR (k — 0.01, right) plotted against the nonlin- 
earity parameter a. Both simulations have been done for 
70 = 10~ 5 , Ao + (u) = 1, and com = 1. The insets show 
the shape of the mechanical potential energy. For negative 
a the springconstant weakens (left), and for a > it stiffens 
(right). 



the oscillator frequency can be dragged along with the 
detuning, thus creating a widely tunable oscillator. 

Finally, it is also possible that the mechanical damping 
of the resonator changes with amplitude as observed in 
the experiments of Rcfs. ||j]-|4l]]. Since in this case only 
the damping rate changes and ojr remains constant, this 
effect can easily be accounted for in the present frame- 
work by sclf-consistcntly solving for the amplitude using 
the average mechanical damping rate. However, as Eq. 
[H] shows, the final amplitude depends only weakly on the 
damping rate of the resonator, so the influence of nonlin- 
ear damping on the maximum amplitude is expected to 
be small. 



VII. CONCLUSIONS 

The backaction limits on optomechanical motion have 
been explored using the amplitude dependent cavity re- 
sponse at the oscillation frequency. Above the thresh- 
old self-sustained oscillations grow until the resonator 
reaches a limit cycle set by the dynamic range of the 
cavity. The largest amplitudes are obtained with an un- 
resolved cavity, but the highest quantum efficiencies are 
found in the resolved case. The latter can be much larger 
than one because many sidebands are involved, result- 
ing in multi-phonon emission. For large amplitudes the 
motion is anharmonic, and numerical simulations show 
that the final amplitude is insensitive to the cavity detun- 
ing and linewidth. Mechanical nonlinearities only have a 
modest effect in the unresolved sideband regime, whereas 
in the resolved case the shift in oscillation frequency due 
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to anharmonicitics can have a strong effect on the ampli- 
tude. 



Acknowledgments 

M.P thanks the Netherlands Organization for Scien- 
tific Research (NWO) / Marie Curie Cofund Action for 
support via a Rubicon fellowship. H.X.T. acknowledges 
support from a Packard Fellowship in Science and Engi- 
neering and a career award from National Science Foun- 
dation. This work was funded by the DARPA/MTO OR- 
CHID program through a grant from AFOSR. We thank 
Florian Marquardt for discussions. 



Appendix A: Numerical simulations 

Eqs. |l| and || are integrated numerically using a fourth 
order Runge-Kutta algorithm with adaptive step size, 



which gives the discretized timetraces of the displace- 
ment, velocity, and (complex) cavity amplitude as out- 
put. The transients where the oscillator relaxes towards 
its final amplitude are discarded so that the steady-state 
oscillations remain. The cavity response <&{A) (Fig. |) 
is found by multiplying the cavity occupation n(t) by 
exp(— iujqt and averaging over many cycles. For the data 
in Fig. U an initial state with a much larger amplitude 
was used to let the oscillator relax preferentially to the 
solution with the highest amplitude among the set of sta- 
ble limit cycles. The amplitude A is obtained by calcu- 
lating the root-mean-squared value of the displacement 
and multiplying by \/2. This ensures that the obtained 
value equals the amplitude for sinusoidal oscillations. An 
alternative definition where half the difference between 
the minimum and maximum value of u{t) is used has 
almost identical numerical values, but naturally shows 
more scatter near the chaotic regions. 
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